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ABSTRACT 

We report Fermi/LAT observations and broad-band spectral modeling of the radio- 
loud active galaxy 4C-j-55.17 (z=0.896), formally classified as a flat-spectrum radio 
quasar. Using 19 months of all-sky survey Fermi/LAT data, we detect a 7-ray con- 
tinuum extending up to an observed energy of 145 GeV, and furthermore we find no 
evidence of 7-ray variability in the source over its observed history. We illustrate the 
implications of these results in two different domains. First, we investigate the origin of 
the steady 7-ray emission, where we re-examine the common classification of 4C +55.17 
as a quasar-hosted blazar and consider instead its possible nature as a young radio 
source. We analyze and compare constraints on the source physical parameters in both 
blazar and young radio source scenarios by means of a detailed multiwavelength analysis 
and theoretical modeling of its broad-band spectrum. Secondly, we show that the 7-ray 
spectrum may be formally extrapolated into the very-high energy (VHE; > 100 GeV) 
range at a flux level detectable by the current generation of ground-based Cherenkov 
telescopes. This enables us to place constraints on models of extragalactic background 
light (EBL) within LAT energies and features the source as a promising candidate for 
VHE studies of the Universe at an unprecedented redshift of z=0.896. 



8 Subject headings: galaxies: active — galaxies: individual (4C +55.17) 

9 — gamma rays: observations — radiation mechanisms: non-thermal 



galaxies: jets 



Introduction 



11 The radio-loud active galaxy 4C +55.17 (0954+556), formally classified as a flat spectrum 

12 radio quasar (FS RQ), has a history of 7-ray observatio ns dating back to the EGRET era, as 

13 3EG J0952+5501 (JHartman et al.lll999l : [Mattox et aklboOlh andEGR J0957+5513 (jCasandiian k Grenier 



20081 ). Due to a relatively poor localization of the EGRET source, however, the association of the 



15 7-ray emitter with 4C +55.17 remained tentative at that time. After the successful launch of the 

16 Fermi Gamma-Ray Space Telescope in June 2008, this associatio n was on the other hand quickly 

17 confirmed by the L arge Area Telescope (LAT: lAtwood et all bood) , initially as OFGL J0957.6+5522 

18 (|Abdo et al.ll2009aH bl^. and most recently as IFGL J0957.7+5523 JAbdo et ahlboiOal ). 



19 The quasar classification of 4C +55.17 may be attributed to the presence of broad optical 

20 emission lines in its spectrum ( Wills et al.lll995l) and high optical/UV core luminosity (absolute 

21 S-band magnitude, Mb < —23; IVeron-Cettv k Veron IboOa). Its redshiflrj , z = 0.896, is based on 



22 the detection of Lya and CIV hues with the HST-FOS (jWills et al.lll995l 1 and Mg II in the SDSS 



i) 



^^ Assuming a ACDM cosmology with Hq = 71 kms ^ Mpc ^, Q.m 
di, = 5785 Mpc, and the conversion scale is 1 mas — 7.8 pc. 



0.27, and Qa = 0.73, the luminosity distance 



23 spectrum (jSchneider et al.ll200?1 ). The optical-UV properties of the source, together with its high 

24 7-ray luminosity of the order L^ ~ lO^^ergs"^, have in turn led to the common classification of 

25 4C +55.17 as a blazar/FSRQ. 

26 However, 4C +55.17 also exhibits a number of morphol ogical and spectral properties that 



27 have placed its exact blazar/FSRQ classification into question (jMarscher et al.ll2002l : iRossetti et al 



12; 



20051 ). FSRQs are uniquely characterized by the presence of a central compact radio core exhibit- 



29 ing a highly variable flat-spectrum co ntinuum, high brightnes s temperatures (Tb), and, typically, 

30 superluminal motions on VLBI scales (jUrrv &: Padovanilll995l ). Indeed, all of the aforementioned 

31 radio properties are shared by the luminous blazars detected in 7-rays: these are exclusively ob- 

32 served to possess compact, highly polarized jets a few milli-arcseconds (mas) in angular size, and 



33 unresolve d radio cores with b rightness tem peratures in the rang e Tb = 10 — 10 K when observed 



114- 



34 at 5 GHz (JTaylor et al.ll2007l ) and 15 GHz (JKovalev et al.ll2009l ). In comparison, 4C +55.17 demon- 

35 strates none of these characteristics. To date, the source shows no evidence of blazar flaring at 

36 any wavelength, nor any evidence of long-term variability, with the exception of a ~ 30% optical 

37 flux-density change noted over a period of 7 years between recent Swift/WOT measurements and 

38 archival SDSS data (see ^ 12.21 and ^ 13.11 for discussion). Furthermore, the VLBI radio morphology 

39 of the source is extende d over ^ 400 pc (p rojected). The peak surface brightness in a VLBA 15 



40 GHz image taken from lRossetti et al.l (J2005l ) is found in the northernmost component and is clearly 

41 resolved, with a corresponding bri ghtness temperature Tb < 2 x 10^ K (consistent with a measure- 

42 ment at 5 GHz; iTavlor et al.l 120071 ). which is uncharacteristic of all the other known quasar-hosted 

43 7-ray blazars. 



44 Based on the radio morphology of 4C +55.17, iRossetti et al.l (J2005l ) firs t suggested t hat the 

45 source may in fact belong to the family of young radio sources (for a review, see lO'Dealll998l ). rather 

46 than blazars. Such sources are characterized by a very low radio variability (if any) and symmetric 

47 double radio structures resembling "classical doubles" on much smaller scales: linear sizes < 1 kpc 

48 for compact symme tric objects (CSOs) and ~ 1 — 15 kpc for medium symmetric objects (MSOs; 



Augusto et al.ll200d ). to be compared with the typical linear sizes of "regular" Fanaroff- Riley type- 
so II radio galaxies of ~ 100 kpc. In many cases, CSO sources are found to exhibit a turnover in their 

51 rad io spectra in the rang e of 0.5— 10 GHz, as the so-called Gigahertz Peaked Spectrum (GPS) objects 

52 do (|de Vries et al.ll 19971 ): similarly, MSO's often displ ay turnover freque ncies below 0.5 GHz, typical 

53 of the Compact Steep Spectrum (CSS) class of sources (JFanti et al.lll99d ). The overlap between CSO 

54 and GPS samples, as well as between sa mples of MSOs and CSS sources, is however not complete 

55 (|Snellen et ahlboool : lAueusto et ahlbood ). In the case of 4C +55.17, the VLBI mo rphology at 5 GHz 

56 reveals two distinct emission regions, to the north and south (JRossetti et al.ll2005l . see also Figured]), 

57 covering a total angular extent of 53 mas (= 413 pc, projected). On the kpc scale, the source reaches 

58 4". 5 (~ 35 kpc, projected), and it is resolved with the VLA in three components, the central one 

59 hosting the VLBI structure. The northern component of the pc-scale emission fea tures a compact 

60 region with a relatively flat spectrum (a = 0.4; F^, oc i/~"; IRossetti et al.ll2005l ). which can be 



61 attributed to a core or a hotspot region, while the southern component features a more diffuse and 



62 slightly steeper-spectrum (a = 0.49) region. iRossetti et al.l (|2005l ) have pointed out that these two 

63 components resemble more compact hotspots and lobes, suggesting a CSO/MSO classification for 

64 this object. The kpc-scale emission might thus be interpreted as a remnant of previous jet activity, 

65 as this is a common feature among sources that show evid ence of intermittent behavior (e.g.. 



Rossetti et al. 



Baum et al.lll990l : ILuo et al.ll2007l : lOrienti k Dallacasall2008l l. Under the CSO/MSO framework. 



found no core candidate between the VLBA-scale lobes at a level > 2 mJy/beam in 



68 a 15 GHz map. 

69 An 11-month comparison of the 7-ray variability and spectral p roperties of 4C +55 .17 against 
the other LAT FSRQs highlights the unusual nature of the source (jAbdo et al.ll2010al b|). Among 



71 all of t he sources originally detected in the 3- month LAT Bright AGN Sample (LB AS; 



Abdo et al 



2009bl'l that were c lassified as FSRQs, 4C +55.17 is characterized by the lowest variability index 



73 (JAbdo et al.ll2010al ). In addition, the unusually hard 7-ray continuum (that is, with a low photon 

74 index F) is found to be one of the hardest among FSRQs in the 1st LAT AGN Catalog (ILAC; 



75 lAbdo et al.M2010bl ). In fact, of those sources included in the ILAC (FSRQ or otherwise) with 

76 > 1 GeV fiux greater than or equal to that of 4C +55.17, only five - all of which are BL Lac objects 

77 (PKS 2155-304, Mkn 421, 3C 66A, PC 1553+113, and PKS 0447-439) - appear with a harder 

78 7-ray spectrum. 

79 In this work, we re-examine the high-energy 7-ray (> 100 MeV) properties of 4C +55.17 using 

80 19 months of LAT all-sky survey data and discuss the implications of these results in two domains. 

81 First, we reconsider the underlying physical processes responsible for the 7-ray emission through 

82 detailed broadband modeling of the source in the context of two scenarios: "young radio source" 

83 and "blazar." In addition, we demonstrate that the unusual properties of the source make it an 

84 ideal candidate for studying the high-redshift universe at very- high energies (VHE) , in particular 

85 for placing constraints on the level of extragalactic background light (EBL). The paper is organized 

86 as follows. Section[2] details the analysis of 19 months of LAT data and discusses the supporting 

87 multiwavelength observations. In particular, section [2?T] focuses on the LAT data reduction, pre- 

88 senting new spatial (localization), spectral, and variability analysis, including a detailed analysis 

89 of the 145 GeV photon detection associated with the source (see also Appendix [A]). Section 12.21 

90 discusses the multiwavelength observations, including analysis of archival radio and Swift X-ray and 

91 optical data, as well as a new hard X-ray detection with the Swift Burst Alert Telescope (BAT). 

92 Spectral properties and classification of 4C +55.17 are discussed in section[3l We follow with a 

93 detailed analysis of the high energy spectrum of 4C +55.17, where we place constraints on models 

94 of EBL and discuss the implications of the 145 GeV photon detection to future VHE observations 

95 of the source ( ^ 13. 2p . Our conclusions are presented in section^ 



96 2. Observations 

97 2.1. FermifLAT Observations 

98 The -Ferm i/LAT is a pair creat ion telescope designed to cover tlie energy range from ~ 20 MeV 



99 to > 300 GeV (jAtwood et al.ll2009l ). The LAT instrument features an improved angular resolution 

100 (^68% = 0-8° at 1 GeV) over previous instruments and a large field-of-view of 2.4 sr. The nominal 

101 mode of operation is an all-sky survey mode, which provides nearly uniform sky coverage approxi- 

102 mately every 3 hours. The following analysis is comprised of 19 months of nominal all-sky survey 

103 data extracted from a 10° region of interest (ROI ) around the J2 000.0 radio position of 4C +55.17 



io4(R.A.= 09''57™38.1844^ Decl. = 55°22^57.769^^: iFev et al.l l2004l ) and covers the mission elapsed 



105 time (MET) 239557417 to 289440000 (August 4, 2008 through March 4, 2010). A 100-second in- 

106 terval at MET 251059717 was removed in order to avoid contamination from GRB 081215A, which 



107 fell within the ROI. Event selections include the "diffuse" event class (lAtwood et al.l |2009| ) rec- 

108 ommended for point source analysis, a zenith angle cut of < 105° to avoid contamination from 

109 the earth limb, and rocking angle cuts at 43° and 52°, respectively, for times corresponding to a 
no change in the instrument's rocking angle from 39° to 50° that occured on September 3, 2009 (MET 

111 273628805). Science Tools v9rl6pl and instrument response functions (IRFs) P6_V3_DIFFUSE were 

112 used for this analysis. 

113 The 19 month LAT localization of 4C +55.17 was determined using gtfindsrc, resulting in 

114 a best-fit position (J2000.0) of R.A. = 09'^57"^40% Decl. = 55°23'40", which is 0.012° = 0.7' offset 

115 from the radio position and falls within the 95% error circle rgg% = 0.017° = 1.0' (statistical 



116 only) . In order to model the 7-ray emission, all point sources from the IFGL catalog (jAbdo et al. 



ii7 l2010al ) within 15° of the source were included. Sources within 10° of the 4C +55.17 radio position 

118 were modeled with their flux and spectral parameters set free, while those sources that fell outside 

119 the 10° ROI were fixed at their catalog values. The diffuse background was modelled using the 

120 recommendeccj Galactic diffuse gll_iem_v02.f it along with the corresponding isotropic spectral 

121 template isotropic_iem_v02.txt. 

122 Prior to fitting the spectrum, the high energy photons attributable to the 4C +55.17 position 

123 (both radio and 7-ray) were found by comparing the energy and incoming angle 9 (defined with 

124 respect to the spacecraft zenith) of each photon within the ROI to the 95% containment radius 

125 of the point spread function defined by the P6_V3_DIFFUSE IRFs. Included among those photons 

126 was a 145 GeV event at an angular separation of 0.06° (R.A. = 09^58^03% Decl. =55°24'00") from 

127 the 4C +55.17 position, falling well within the 95% containment radius for the given energy and 

128 angle of incidence. Through an analysis of the event diagnostics, the photon nature of this event 

129 is confirmed here for the first time (for further details regarding the 145 GeV event analysis, see 

130 Appendix [A|). In addition, several photons in the ~ 30 — 55 GeV range were also detected. The 

131 association of the 145 GeV photon with 4C +55.17 tentatively places it as the highest-redshift source 



']http://feriiii.gsfc.iiasa.gov/ssc/data/access/lat/BackgroundModels.html| 



32 to be observed at VHE to date. 

33 A spectral analysis of 4C +55.17 was performed with gtlike using the LAT data between 

34 100 MeV and 300 GeV. Spectral data points were first obtained by fitting each of 9 equal logarith- 
ms mically spaced energy bins to a separate power law with index and prefactor parameters set free. 

36 From the resulting data points, a break in the spectrum could be seen to occur at ~ 1.6 GeV. This 

37 was confirmed by performing an independent unbinned likelihood fit over all the data from 100 MeV 

38 to 20 GeV using power law, log parabola, and broken power law models, with the break energy of 

39 the broken power law fixed at the peak in the vF^, representation (Sbr ~ 1.6 GeV). The maximum 

40 energy of 20 GeV was chosen in order to avoid fitting any portion of the spec trum that may be 

41 significantly attenuated by the EBL. A likelihood ratio test (JMattox et al.lll996l ) resulted in a 4. la 

42 improv ement of the broken powe r law over the single power law used in previous analyses of the 

43 source (lAbdo et al.ll2009al . l2010al ). as compared to a 3.8a improvement over the power law from the 

44 log parabola. We therefore consider the broken power law to be the most accurate representation 

45 of the intrinsic 7-ray spectrum of the source. 

To test the 7-ray variability over the 19-month period, we made light curves in time bins of 7 

47 and 28 days. Due to the limited statistics over each interval, the source was fit to a single power- 

48 law in each bin, with index and prefactor parameters free. To improve the fit conv ergence, point 



49 sources in the ROI were included only if they were detected with a test statistic (TS: lMattox et al. 
19961 ) greater than 1 (~ lo"). The resulting light curve (> 100 MeV), divided into 7 day bins. 



51 is shown in Figure[2j The variability of 4C +55.17 was analyzed by means of a x^ test, where 

52 we assumed the model describing the data to be a constant straight line with intercept equal 

53 to the weighted mean of all > 3a detections. This test yielded a x^ probability P{x^ > Xobs) 

54 of 0.96 and 0.87 for the 7 day and 28 day light curves, respectively, and was thus in agreement 

55 with the tested hypothesis. We therefore found no evidence of variability in 7-rays over the 19- 

56 month LA T observing period, c onsistent with the previous 11-month lightcurve analysis (~ 30 day 

57 bins) from lAbdo et al.l (J2010al ). In addition, the weighted mean for this period was found to be 



58 (9.5 ± 0.4stat + 0.83sys — 0.49sys) X 10 ^ ph cm ^ s ^ , which is consistent with the EGRET measured 



59 flux of (9.1 ± 1.6) X 10^^ phcm^^s^^ ( Hartman et al.lll999l ) as well. Systematic uncertainties on 

60 the LAT flux were determined by bracketing the instrument effective area to values of 10%, 5%, 
.61 and 20% their nominal val ues at log (E /M e V) = 2,2.75, and 4, respectively. We note that these 



62 flndings differ from those of lNeronov et al.l (I2OIOI ). who claim variability between the EGRET and 

63 LAT measured fluxes. We believe this discrepancy lies in a mis-quoted value of the EGRET flux. 



2.2. Multiwavelength Data 

2.2.1. X-ray 



166 We analyzed all Swift (JGehrels et al.ll2004l ) data obtained over the 19 -month LAT observing 

167 period, which consisted of three X-ray Telescope (XRT: lBurrows et al.ll2005l ) snapshots (1.6-4.5 ks). 



168 in order to check the X-ray state of the source. We used the xrtgrblc script (available in the HEASoft 

169 package version 6.8) to analyze the XRT observations: we reprocessed the data stored in the 

170 HEASARC archive using the latest XRT calibration database (20091130), selecting the events with 

171 0-12 grades in photon counting mode (PC). The scripts chose the optimal source and background 

172 extraction regions based on the source intensity: the X-ray photons were extracted using a 25" 

173 circle for the source and an annulus with 50" — 150" inner-outer radius for the background. Adding 

174 all of the exposure and performing a C-statistic fit from 0.3 — lOkeV using XSpecl2, we found 

175 the best fit obtained to be a power law with absorption fixed at the galactic value (A^h = 9 x 



176 10 cm ), where we obtained the photon index T = 1.84ib0.19, with an absorbed flux of (8.3 



-ia) 



10 ^^ erg cm ^ s ^ and an unabsorbed flux of 



^.5+}:I)x 10-13 



erg cm ^ s ^ . 



Comparing each of the 

178 individual observations, no X-ray variability was found, with all measurem ents falling within the 

179 joint errors. These results were also compared with previous Chandra data (jTavecchio et al.l 120071 ) 

180 obtained June 16, 2004, where the flux was found a gain to be non-variab le within the statistical 

181 errors. Finally, historical X-ray data from ROSAT (jComastri et al.l 119971 ) obtained November 7, 

182 1993 were included in the spectral energy distribution (SED) modeling to further constrain the soft 

183 X-ray portion of the spectrum. 



184 In the hard X-rays, data from the Swift/HAT ([Aiello et al. 



2008 



20091 ) were analyzed using 



185 five years of cumulative exposure from November 2005 — 2010. We detect the source for the first 

186 time in the hard X-ray band, with a 15 — 150 keV flux of (6.75l^5'2i) x lO^^^ ergcm"^ s~i and a 



187 power-law photon index, T = 1.79 



+ 1.17 
-0.84- 



2.2.2. Optical & Infrared 



189 During each o f the three Swift pointings in 2009, Ultra- Violet/Optical Telescope (UVOT; 



Roming et al.l l2005l ) observations were also obtained. Data were obtained in all 6 filters in the 



191 first two epochs, and the last epoch with only the W2 filter. The data reduction and analysis was 

192 performed using the uvotgrblc script, which reprocesses the data stored in the HEASARC using 

193 the latest UVOT calibration database (20100129). The optimal source and background extraction 

194 regions were a 5" circle and a 27"— 35" annulus, respectively. Tabled] summarizes these observations. 

195 A comparison of the results between each epoch shows the source to fall within the joint errors in 

196 flux in the optical to UV bands across all three epochs. These results we re also compared with 

197 archival SDSS data from February 2, 2002 (jAdelman- McCarthy et al.l 120081 ). A comparison of the 

198 UVOT and SDSS [/-band flux densities shows an increase from (0.187 ± 0.003) mJy in the SDSS 

199 data to (0.250 it 0.007) mJy in the UVOT data, indicating a ~ 30% rise in flux over 7 years. In 

200 addition, the UVOT V- and B-hand flux densities were averaged using a least-squares approach to 

201 a linear flt and compared with the SDSS 5r-band, which fell between the two. The average of the 

202 UVOT V- and i?-bands, measured at (0.305 it 0.014) mJy, shows a similar ~ 25% increase from 

203 the SDSS measured value of (0.240 it 0.011) mJy. A comparison of the Swift UV OT measurements 

204 to the continuum flux underlying the Lya line obtained by HST-FOS in 1993 (jWills et al.lll995l ) 



205 shows the fluxes to be equal between these two periods. 



206 In the near-infrared, we included historical data from the 2MASS Point S ource Catalog (|Cutri et al. 



20031 ). for which the absolute calibration was taken from lCohen et al.l (J2003l ) . All infrared, optica l, 



208 and ultraviolet data were dereddened by means of the extinction laws giv en bvlCardelli et al.l (119891 ) 



209 assuming a i?-band Galactic extinction (^b = 0.038) as determi ned vialSchleeel et al.l (119981 ). and 



210 a ratio of total to selective absorption at V equal to Ry = 3.09 (jRieke &: Lebofskvlll985l ). 



2.2.3. Radio 



212 To model the 7-ray emi ssion in 4C+55.17 (sec. 3.1), we compiled integrated ra.dio to sub-mm 



213 measurements of the source (Bloom et al 



1994 



Huang et al.lll998l : iReich et al.lll998l : Ijenness et al 



2i4 l20ld ). including 5-year WMAP data (jWright et al.ll2009l ). and other archival data from the NASA/IPAC 
215 Extragalactic Database (NED). In order to isolate the total radio flux from the inner ~ 400 pc 



216 scale structurCj, we re-analyzed several archival VLA data sets from 5 to 43 GHz (see Figures [3] 

217 andSj). The typical resolutions are ~ 0.1" to 0.4", ensuring a total me asurement of th e ~ 50 mas 



218 scale structure without loss of flux as in the VLBI observations (e.g., iRossetti et al.l l2005l) . We 



219 also include similar measurements from p reviously published VLA 5 and 8. 4 GHz (IReid et al. 



Mvers et al. 



2003 



Tavecchio et al.ll2007l ) and MERLIN 0.4 and 1.7 GHz (JReid et al.lll995l ) maps. 



1995 



221 The radio variability properties of 4C +55.17 are important for assessing its nature. We 

222 therefore search ed the literature for various archival radio to sub-mm monitoring observations of 



223 the s ource (e.g. lAltschuler &: Wardldll976l : IWardle et al. 



1981 



Seielstad et al. 



1983 



20 id ) , including 22 and 37 GHz data from the Metsahovi monitoring program ( Terasranta et al. 



Jenness et al. 



199 



20041 . I2OO5I ) . While the I War die et al.l (|l98ll ) data was not directly available, we note from the 



226 literature that the authors found the source to be non-variable. Variability in each of the remain- 

227 ing cases was measured by applying a statistical x^ test of the available data using the hypothesis 

228 of a constant source with flux equal to the weighted mean. The results were consistent with the 

229 tested hypothesis in each case, with the exception of the Metsahovi data, which yielded probabilities 

230 P(x^ > Xobs) of 6.44 X 10"^^ and 8.56 x IQ-^"^ at 22 GHz and 37 GHz, respectively. To quantify this 

231 variability, we compared fractional variabilit y indices using the formu la VarAs = (Smax — Smin)/Smin 

232 used in a variability study of GPS sources (JTorniainen et al.l 120071 ) . where we obtained values of 

233 3.5 and 1.43 at 22 and 37 GHz, res pectively. The 2 2 GHz value fell slightly above the nomi- 

234 nal variability threshold of 3.0 set by iTorniainen et al.l (|2007l ) as an upper limit for the bona fide 

235 GPS sources. This result, however, arose due to a single outlying flux measurement at 22 GHz of 

236 0.32 lb 0.09 Jy which occurred ~ 40 minutes after a previous measurement of 1.12 it 0.08 Jy at the 



^^The kpc-scale radio emission is not expected to contribute significantly toward the modeling of the high energy 
portion of the spectrum (see JJS.l.ll fc ^3.1.2[) . 



9 



237 same frequenc}0. Removing this questionable flux point and performing the test again resulted 

238 in a fractional variability index of 0.89, which fell well within the proposed threshold for genuine 

239 GPS galaxies. We therefore find the degree of variability in 4C +55.17 to be consistent with the 

240 behavior of confirmed young radio galaxies, rather than blazars. 



Results 



242 3.1. Modeling & Classification 

243 3.1.1. CSO Modeling 

244 As noted in the introduction, there are several reasons to consider the possible nature of 

245 4C +55.17 as an example of a luminous AGN exhibiting recurrent jet activity, with young and 

246 symmetric (CSO-type) inner radio structure instead of a "core-jet" morphology typical of blazars. 

247 While the physical nature and the origin of the CSOs is at some level still debated, the most likely 

248 and widely a ccepted hypothesis states that they are the young versions of present-day extended ra- 



249 dio galaxies (jPhilips fc Mute] 



1982 



Fanti et al.lll995l ). In the alternative explanation, these sources 



250 are considered to be of a similar age to normal radio galaxies, but only co nfined /frustrated due to 

251 dramatic interactions with a surrounding dense gas in their host galaxies (Ivan Breugel et al.lll984l : 



Wilkinson et al.l Il994l ) . The latter scenario is however inconsistent with the lack of observational 



253 evidence for th e amount of ambient gas required to supply suffic ient confinement (|De Youndll993l : 



Carvalhd Il994l . 119981 : ISiemiginowska et al 



20051 : see, however, iGarcia-Burillo et al 



20071 for no- 



255 table exceptions). More promising is the refore the "youth" scenario for C S Os, for which a number 

256 of evolutionary models were proposed (jBegelmanl Il996l : iDe Youna 119971 : iPerucho &: Martil |2002| : 



Kawakatu fc Kindl2006l ) 



258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 



While many observational properties of 4C +55.17 make its classification as a young radio 
source compelling, it is also worth noting the characteristics that could make such a classification 
potentially difficult. For example, if 4C +55.17 is indeed a CSO, it is the only such object to be iden- 
tified as a 7-ray emitter in IFGL/ILAC, with a GeV fiux nearly an order of ma gnitude higher tha n 
the lower limit of the complete fiux-limited subsample within the ILAC catalog (|Abdo et al.ll2010bl ). 
This would immediately set the object apart as an outstanding member of its class. In addition, the 
relatively high radi o polarization of the source (~ 3% in a ~ 0.2" resolution VLA 8.4 GHz image; 



Jackson et al.ll2007ll. is unch aracteristic of the typically low (< 1%) radio polarization seen among 



CSOs (IReadhead et al.lll996l ). although polarized emission from CSOs has occasionally been found 



(e.g., iGugliucci et al.l 120071 ). The low polarization of CSOs, which are entirely embedded within 



the inner regions of the host galaxy, is often attributed to the large expected Faraday depths of 



^^Variability within hour timescales is rare at the frequencies observed by Metsahovi (A. Lahteenmaki, T. Hovatta, 
& M. Tornikoski, private communication 2010) 
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269 the surrounding interstellar medium (JBurn 1119661 : iBicknell et al.lll99?1 : lGugliucci et al.ll2007l ). The 

270 surrounding medium may also play a key role in shaping the spectral turnover seen in the GPS class 

271 of young radio sources, through the free-fre e absorption (FFA) process (either internal or external 

272 to the emission region; IBicknell et alj 119971 : iBegelmanI Il999l : iPeck et al.lll999l ). The nature of the 

273 absorber is however still widely debated, ar id both the synchrotron self-absorption (S SA) and FFA 

274 processes are considered as viable options (jO'Dea &: Bauinlll997l : ISnellen et alJl200Cl ). 



275 If FFA effects are indeed responsible for the spectral turnover in GPS sources, then the rel- 

276 atively flat (a ~ 0.4 — 0.5) power-law radio continuum of 4C -1-55.17, which shows no indication 

277 of a low-energy turnover, may indicate an exceptionally small amount of ionized ambient gas in 

278 the vicinity of its young radio structure. More specifically, if the radio absorber may be identified 

279 with ionization-bounded hydrogen clouds of interstellar matter pr esent at pc to kp c distances from 

280 th e center and engulfed by the expanding lobes, as proposed by iBegelmanI (119991 ) and advocated 

281 bv IStawarz et al.l (J2008l ). and if a significant part of this gas has been evacuated prior the onset 

282 of new jet activity, then one would expect much less severe absorption of the low frequency radio 

283 emission, resulting in a lower turnover frequency compared to that of GPS galaxies. In this case, 

284 the relatively high polarization of 4C -1-55.17 (as for a young radio source) would find a natural and 

285 straightforward explanation as well. 

286 In considering the hypothesis outlined above, and in order to investigate the 7-ray emission 

287 detected from 4C -1-55.17 in a framework that is more consistent with the observed properties 

288 of the source, we ap ply the dynamical model for the broad-band emission of CSOs proposed by 



Stawarz et al.l (l200a) and successfu l ly test ed against a sample of X-ray detected young radio galaxies 



290 of the CSO type bv lOstorero et al.l (J2010l ). In this model, the newly born relativistic jets propagate 

291 across the inner region of the host galaxy and inject ultrarelativistic electrons into the compact 

292 lobes. These electrons, which provide the bulk of the internal lobes' pressure, cool radiatively and 

293 adiabatically within the sub-relativisticly expanding plasma, thus producing isotropic synchrotron 

294 (radio) and IC (X-ray to 7-ray) radiation. In the model, the broad-band emission spectra are 

295 evaluated self-consistently for a given set of the initial parameters of the central engine and of the 

296 host galaxy, taking into account the time-dependent evolution of the radiating electrons. For a 

297 given linear size of the system, which is uniquely related to a particular age of the system, the 

298 observed broad-band emission spectrum is given as a snapshot of the evolving multiwavelength 

299 radiation of the lobes. Based on this model. IStawarz et al.l (120081 ) argued that, in fact, young radio 

300 galaxies should be detected by Fermi/LAT at GeV photon energies, albeit at low fiux levels and 

301 after an exposure longer than one year. Other (physically distinct) scenarios for the production 

302 of soft, high energy, and VH E 7-rays in the lobes an d hot spots of young r adio galaxies have been 

303 proposed and investigated bv iKino et al.l (|2007l . l2009l ) and lKino fc Asand (|2010l ). 



304 In the more detailed description of the model, the jets with total kinetic power (Lj) propagate 

305 with the advance velocity (vh) in the interstellar medium, characterized by a given number density 

306 (n-ext)- At a particular instant of the source evolution, the infiated lobes will have a corresponding 

307 linear size {LS). The electrons injected through the termination shock into the lobes with the 
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308 intrinsically broken power-law energy distribution cool due to the synchrotron and IC processes. 

309 The most relevant ambient photon fields for the IC scattering are the UV emission of the accretion 

310 disk (mean photon energy Sdisk = 10 eV, disk luminosity Ldisk), the starlight (ffstar = 0.83 eV, 

311 host luminosity Lstar), and the infrared emission of the obscuring nuclear torus (ffdust = 0.02 eV, 

312 dust luminosity L^^gi). The magnetic field intensity is expressed in terms of the ratio of energy 

313 densities stored in the radiating electrons and the magnetic field, U(./U-q, which is constant during 

314 the expansion of the radiating plasma. Note, however, that [4 and U-q, as well as the energy 

315 densities of the ambient photon fields (and hence the elect ron cooling conditions) do change with 



316 time, and therefore depend on LS (see lStawarz et al.ll2008l . for more details). 



317 The fit of the "young radio source" model to the collected broad-band dataset for 4C -1-55.17 

318 is illustrated in FigureO In fitting the SED, we assume that the projected source size of the inner 

319 radio structure {LS ~ 400 pc) is equal to the actual source size, though we note that this may be 

320 underestimated due to possible projection effects. Indeed, some amount of projection off the plane of 

321 the sky is required to account for the presence of the intense disk-related optical/UV continuum and 

322 the broad optical emission lines in the spectrum of 4C +55.17 (which might otherwise be completely 

323 obscured), as well as to account for the asymmetry in brightness between the two lobes. In fitting 

324 the broadband SED, the following model-free parameters were obtained: Lj ~ 6.6 x lO^^ergs"^, 
325Ldisk ~ 2 X 10^6 erg s-i, Lstar ^ 10^5 erg s-i, Ldust - 10^5 erg s'^, UJUb ~ 160, v^, ~ 0.3c, and 

326 next — 0.1 cm^'^. The injection electron energy distribution is characterized by the minimum, break, 

327 and maximum electron Lorentz factors, 7min — 1, 7br — 2 x 10^, and 7max — 4 x 10^, respectively, as 

328 well as by the low- and high-energy electron spectral indices, si ~ 0.5 and S2 — 2.5. The model fits 

329 quite successfully all the relevant data points within the low-frequency (radio) and high-frequency 

330 (hard X-ray to 7-ray IC component) ranges; it also reproduces nicely the spectral break within the 

331 Fermz/LAT photon energy range. We note that in our modeling here and below we do not consider 

332 7-ray absorption effects related to the direct or reprocessed emission of the accretion disk, which 

333 may lead to the attenuation of the lobes' (or jets') emission at photon energies > 100 GeV (see in 



334 this context Tanaka et al. 2011 



335 Looking closely at the UV part of the spectrum, we note an approximate factor of two differ- 

336 ence between what is observed and what is required for producing the appropriate luminosity in 

337 IC-scattered 7 rays. This can be resolved by recalling that in the framework of the model the opti- 

338 cal/UV photon energy range is dominated by the thermal UV disk emission that may suffer from 

339 some non-neglible obscuration by the circumnuclear dust for moderate inclinations of the source to 

340 the line of sight. Also worth noting are the variation timescales of the di sk, which are governed b y 

341 the viscous motion within tens of gravitational radii from the black hole ( Collier &: Petersonll200ll ). 



342 This can account for the ~ 30% variation over seven years seen between the optical measurements 

343 from UVOT and SDSS (see ^ 12. 2p . On the other hand, the CSO-related non-thermal IC emission is 

344 expected to be non-variable in accordance with the observations, because this emission is produced 

345 within the hundred-pc-scale and sub-relativistictically expanding lobes, and hence the UV photons 

346 seen by the lobes' electrons will be averaged over the entire spatial extent of the radio structure. 
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347 Here we do not model the accretion-related emission in detail, but only roughly represent it as a 

348 blackbody component for the purpose of the evaluation of the IC radiation of the lobes. Likewise, 

349 the steep-spectrum soft X-ray continuum is not accounted for by the IC emission of compact lobes 

350 and instead may be attributed to the radiativ e output of the accretion disk and its corona (see 



Siemiginowska et al.l 120081 : ISiemiginowskal l2009l . for the X-ray properties of young radio sources) 



352 Yet it should be also noted that the particular CSO model presented here cannot account for the 

353 millimeter-to-near infrared emission of 4C +55.17. In the framework of the discussed scenario, this 

354 has to be attributed to the radiation of the underlying jet, and not of the compact lobes. 

355 The physical parameters of 4C +55.17 emerging from the model fit presented above may be 

356 compared with the ph ysical paramete r s of b ona fide young radio galaxies derived in the framework 

357 of the same model by lOstorero et al.l ()2010l ). The most significant differences can be noted in the 

358 kinetic luminosity of the jet (Lj), the UV luminosity of the accretion disk (Ldisk), and the electron- 

359 to-magnetic field energy density ratio {Ue/U-o). In particular, the jet and the disk luminosities 

360 of 4C +55.17 are higher (by one to two orders of magnitude, on average) than the analogous 

361 luminosities of GPS radio galaxies. This is in fact expected, since t he analyzed source is much more 

362 powerful than the relatively low-power radio galaxies modeled by lOstorero et al.l (|20ld ) . The disk 

363 luminosity obtained from the fit can also be compared with t he expected value ba sed on the total 

364 luminosity of emission in broad lin es (Lrt,r). Using e g. (1) in lCelotti et al.l (119971'). along with th e 

365 line fiuxes of 4C +55.17 obtained in I Wills et al.l (|l995l ) and the line ratios from lFrancis et al.l (|l99ll ). 

366 we estimate the value of Lblr to be 1.2 x 10^^ ergs~^. Using the approximation Ldisk — 10 x Lblr, 

367 we thus obtain -Ldisk — 1-2 x lO^^ergs^^, which again falls within a factor of two of the value 

368 obtained through the model, consistent with the level of uncertainty expected using this method. 



369 3.1.2. Blazar Modeling 

370 As already noted in the introduction, the lack of pronounced variability and resolved VLBI 

371 structure in 4C +55.17 would make it a highly unusual case of a blazar/FSRQ. Still, it is a worth- 

372 while exercise to consider the physical parameters implied from the blazar model. In the framework 

373 of the blazar scenario the observed non-thermal emission of this source, including the 7-ray fiux 

374 detected by Fermi/hAT , is expected to o riginate in the inn ermost parts of a relativistic jet that is 

375 closely aligned with the line of sight (e.g., ISikora et al.lll994l ). In this case, the broad-band emission 

376 of 4C +55.17 should be strongly Doppler boosted in the observer rest frame, and variable on short 

377 (days to weeks) timescales. The expected size of the blazar emission region (sub-pc), which is orders 

378 of magnitude smaller than the linear size of the resolved inner radio structure discussed previously 

379 (~ 400 pc), as well as the presence of relativistic beaming effects, constitute the main differences 

380 between the "blazar" and "young radio source" scenarios. 

381 In order to model the broad-band spectrum of 4C +55.17 as a blazar emission, we apply the dy 



382 namic al model BLAZAR developed by lModerski et al.l (120031 ) and later updated by lModerski et al. 



383 (J2OO5I ) for the correct treatment of the Klein-Nishina regime (for applications of the model, see e.g. 
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384 ISikora et al.l l2008l : iKataoka et al.1 12008| ) . The model describes the production of the non-thermal 

385 emission by ultrarelativistic electrons, which are accelerated in situ within thin shells of plasma 

386 propagating along a conical relativistic jet (bulk Lorentz factor, Fj ^ 1, jet opening angle 9j ~ 1/rj) 

387 and which carry a fraction Lg/Lj of the jet kinetic power. The acceleration process is attributed to 

388 the Fermi mechanism operating at strong shocks that are formed within the outflow as a result of 

389 the shells' collisions, which take place at distances greater than rg from the jet base, resulting in the 

390 injection of a broken power-law electron energy distribution into an emission region of linear size 

391 R and magnetic field intensity B. The non-thermal emission evaluated at r ~ -R/^j ^ '''o includes 

392 the synchrotron and IC components, with the target photons for the inverse-Compton scattering 

393 provided by the jet synchrotron radiation and the external photon fields (predominantly accretion 

394 disk emission reprocessed in the broad line region and within the dusty torus) . 

395 The BLAZAR fit to the broad-band spectrum of 4C +55.17 is shown in Figure HI The fit was 

396 obtained with the following free parameters of the model: Lj > L^ ~ Q x 10"^^ erg s^^, i^disk — 

397 3 X lO''^ ergs"\ Ldust =^ 6 x lO''^ ergs"\ ro ~ 4 x 10^^ cm, r ~ 8 x 10^*^ cm, Fj ~ 12, and B ~ 0.2 G. 

398 For the injection electron energy distribution, the electron Lorentz factors 7min — 1, 7br — 1-5 x 10'^, 

399 and 7max — 10^ were obtained, along with the spectral indices, si ~ 0.5 and S2 — 2.8. The blazar 

400 model fit to the collected dataset, and the implied physical parameters of the 4C +55.17 jet and 

401 its central engine, should be regarded as plausible. Notable differences with respect to the CSO 

402 model discussed previously can be however noted within the radio-to-X-ray frequency range. In 

403 particular, unlike the CSO fit, the blazar model fit does not account for the bulk of the observed 

404 radio fluxes. These emissions, in the framework of the blazar scenario, must therefore be produced 

405 further down the jet, at relatively large distances from the blazar emission zone. On the other 

406 hand, the high-energy tail of the synchrotron blazar emission dominates the radiative output of 

407 the system around the observed near-infrared and optical frequencies, and also at soft X-rays. The 
403 observed hard X-ray spectrum of 4C +55.17 can be hardly attributed to the IC blazar emission 

409 and requires an additional spectral component. In general, the CSO and blazar fits differ the most 

410 within the near infrared and X-ray domains, hence future constraints on the hard X-ray and near 

411 infrared spectra, along with continued monitoring from the radio to the 7-ray band, should be 

412 considered as a potential way of discriminating between the two scenarios. 

413 In comparing these two models, we also note the important difference between the blazar 

414 and CSO model for 4C +55.17 in the radiative efficiency of the emission zone. Compact emission 

415 zones of blazar sources are typical ly characterized by a very low (less than a few percent) radiative 



416 efficiency (e.g., ISikora et al.lll994l ). In this context, only a small fraction of the jet kinetic power 

417 is dissipated in the blazar emission zone and radiated away in the form of high-energy emission, 

418 which is strongly Doppler-boosted in the observer frame due to the relativistic bulk velocity of the 

419 emitting plasma. This is also the case for 4C +55.17 when modeled in the framework of the blazar 

420 scenario discussed above. On the other hand, the radiative efficiency of the su b-relativistically 



421 expanding lobes of young radio sources is known to be large, often exceeding 10% (JDe Young 



1993 



Stawarz et al.ll2008l ) , which naturally accounts for the particularly high intrinsic radio luminosity of 
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423 these sources, being comparable to the most powerful radio galaxies and quasars (JReadhead et al 



19961 ). Likewise, when modeling 4C +55.17 as a CSO, the radiative efficiency was similarly high. 



425 The improved radiative efficiency of CSO sources, together with the relatively high jet kinetic power 

426 implied by the young radio source scenario (higher than that implied by the blazar model), can 

427 thus account for the observed 7-ray luminosity even in the absence of relativistic beaming. 

428 While the CSO-type and blazar modelings of the broad-band spectrum of 4C +55.17 can both 

429 account for the 7-ray emission from the source, we find the implied value for the bulk Lorentz 

430 factor Fj ~ 12 under the blazar scenario difficult to reconcile with its observed VLBI properties. 

431 The physical mechanism responsible for the steady 7-ray emission is also not easily explained 

432 under this framework. Still, the unusual characteristics of 4C +55.17 as for a young radio source 

433 may be evidence for a combination of radiation produced in the sub-pc scale relativistic jet and 

434 the emission of the compact lobes. The modeling of this complex scenario, which might require a 

435 combination of the two models discussed above, is beyond the scope of the present work. A similar 

436 situation was recently considered by Migliori (in prep), who have studied the high-energy (X-ray to 

437 7-ray) emission of radio-loud quasars with CSO-type inner radio morphology, such as, e.g., 3C 186. 

438 Objects of that type might be very common in scenarios of intermittent j et production in active 



439 galax i es, proposed to account for the evolution of r adio- loud AGNs (e.g., [Reynolds &: Begelman 



19971 : ISiemiginowska et al.l 120071 : ICzernv et al.l |2009| . and references therein) . With its complex 



441 radio structure featu ring inner and outer lobes, as well as jet-like features (JRossetti et al. 



2005 



Tavecchio et al.l 120071 ) . 4C +55.17 might thus be another example of AGN with intermittent jet 



443 production. 



444 3.2. High Energy 7-ray Continuum of 4C +55.17 

445 At energies > 10 GeV the 7-ray continua of high-redshift sources begin to suffer from sub- 

446 stantial attenuati on by the still poorly known EBL photon field due to the photon-photon pair 



1 poorly 

447 creation process (JHauser &: Dwekl I2OOII ) . By attributing the attenuation of AGN 7-ray spectra to 

448 these interactions, it is thus possible to p lace significant upper limits to the EBL provided some 

449 estimate of the source's intrinsic spectrum ([Aharonian et al.ll2006l ). In this respect, combined Fermi 

450 and VHE measurements by Cherenkov telescopes such as MAGIC, H.E.S.S., and VERITAS, con- 

451 tinue to prove succes sful at providing these limits (e.g. iGeorganopoulos et al.ll2010l : lAleksic at al. 



2011 



Orr et al. 



2011 



Furthermore, with the VHE detection of the FSRQ 3C279 {z = 0.536) by 



453 MAGIC (lAlbert et al.ll2008l ). and the recently announ ced detections of others quasars - PKS 1510- 

454 08 9 (z = 0.361) by H.E.S.S. dWagner fc Beheralboiol ;) and PKS 1222+216 (z = 0.432) by MAGIC 



455 ( Aleksic at al. 



2011 



the search for increasingly distant luminous sources in the observable range of 

456 ground-based Cherenkov Telescopes has become one of considerable interest to the TeV community. 

457 The extension of the observed 7-ray spectrum of 4C +55.17 up to energies of 145 GeV, coupled 

458 with the source's relatively high redshift of 2; = 0.896, immediately places it among the most 

459 important high- 2; objects that can be used for constraining the widely debated EBL level even within 
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460 LA T energies; f o r an ov erview of different methods for constraining the EBL with the Fermi/LAT, 



461 see 



Abdo et alj (J2010d ). Figure[5] illustrates the r^^ opacity at the redshift z = 0.896 due to 7-ray 



462 absorption with the EBL intensity and spect r al distribution for various models (IFinke et alJ 120101 : 



Franceschini et alJ 120081 : iGilmore et al.1 l2009l : iKneiske et alJ |200J; IStecker et al.l 120061 ) considered 



J 



464 as a function of photon energy. The highest-energy photon associated with 4C +55.17 is also 

465 indicated. As illustrated in the figure, attenuation due to the EBL-related absorption of 7-rays 

466 within the observed range is pred icted in all the scenarios , including those c l ose to the lower limi ts 

467 derived from galaxy counts (e.g., iFranceschini et al.ll2008l : iFinke et al.ll2010l : iGilmore et al.ll2009l ). 



468 To test the validity of particular models o f the EBL us i ng the 4C +55.17 spectrum, we followed 

469 the likelihood ratio test method described in lAbdo et al.l (|2010d ). The full > 100 MeV observed 

470 spectrum was firs t fit to a broken power law with EBL attenuation from 9 separat e EBL models 



471 (Finke et al. 



20061 : 



201c : 



199; 



Franceschini et al. 
Salamon h Steckeri 



2008: 



Gi 



Kneiske et al. 



more et al. 



20091 : 



Primack et al. 



200,51 : 



Stecker et al. 



2004 ) . with the normalization of the attenuation 



473 parameter t^^{E,z = 0.896) fixed to 1 at all energies. The results from each of the spectral fits, 

474 including the low (Fi) and high (r2) broken power law indices, as well as the integral flux values, are 

475 summarized in Table[2j Allowing the normalization of the predicted opacity r-y^ to remain free, we 

476 then compared each result with the likelihood values obtained when the normalization parameter 

477 was fixed to 1. In cases where the r^^ normalization was reduced, a rejection at the level of n 

478 standard deviations (a) of the particular model could be established using the formula: 



n 



V-2 X [log (L 



fixed J 



log (Lfree)] , 



(1) 



479 where Lfixcd and Lfree are the likelihood values of the fits for fixed and free normalizations on 

480 T^^, respectively. Using thes e results , we w ere able to reject two separate models at > 3a signif- 

481 icance. These were the iStecker et al.l (J2006l ) baseline and fast evolution models at 3.9a and 4.3(t, 

482 respectively, with preferred ii o rmaliz ations of 0.17 + 0.14 and 0.16 + 0.12. These models were sim- 

483 ilarly rejected in I Abdo et al.l ( 2010a ) by applying the likelihood ratio test to several blazars and 

484 gamma-ray bursts with redshifts rangi ng from z = 1 . 05 to z = 4.24. Combining this result with 



485 th e overall rejection s ignificance of the IStecker et al.l (|2006l ) baseline model of 11. 4a as calculated 



Abdo et al.l (|2010d . §3.2.3 therein), we obtain a new combined rejection of 11.7a for both the 



486 m 

487 baseline and fast evolution models. 

438 Figured] shows the predicted shape of the int rinsic spectruin of 4 C +55.17 obtained by de- 

489 absorbing the observed Fermi spectrum using the IStecker et al.l (J2006l ) baseline EBL model. A 

490 common feature occurring from models which over-predict the level of EBL is that of an unbounded 

491 exponential spectral rise at highest energies - a behavior which can largely be considered non- 
492 physical, and has thus been used in previou s studies to place constraints on the EBL using TeV 
493 obs ervations (e.g . 



494 the Stecker et al 



Dwek &: KrennrichI |2005| ) . This behavior is clearly illustrated in the case of 
(I2OO6I ) baseline model. Such a feature would in turn require the modeling of an 



495 additional spectral component beyond that which we consider in ^ 13.11 and that would be orders of 
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496 magnitude more luminous than the observed inverse-Compton (IC) peak. We also note that any 

497 intrinsic absorption that may be taking place within the source represents an even greater rejection 

498 of t his model, as the tr ue attenuation due to the EBL would be even less. We therefore consider 

499 the lStecker et al.l ( 20061 ) baseline and fast evolution-^ 1 models to over-predict the true level of EBL 



500 at the observed redshift and energies. 

501 With its excellent sensitivity in the high-energy range, the LAT instrument provides a unique 

502 opportunity to search for VHE candidates at high redshifts through detailed spectral analysis of the 

503 Fermi data. In the case of 4C +55.17, the attenuated high-energy spectrum obtained from fitting the 

504 nine often discussed EBL models is illustrated in Figure [71 Each spectrum is extrapolated beyond 

505 the highest observed photon energy of 145 GeV and compared against the typical differential flux 

506 sensitivity curves of currently operating TeV telescopes. With the exception of the four "highest- 

507 level" EBL models (including the two models ruled out by the present work), the observed 4C +55.17 

508 spectrum is found to lie at the observable threshold for ground-based observations. It is also worth 

509 noting that while intrinsic absorption from interactions with the UV disk and infrared torus may 

510 contribute to the spectral attenuation at energies > 100 GeV, this effect would be reduced in cases 

511 where the 7-ray emission takes place at hundreds-of-parsecs scale distances from the central black 

512 hole, for which there is compelling evidence in the case of 4C +55.17 (see §[3TLTJ. In addition, 

513 with the present analysis we find no evidence of variability in 4C +55.17 over 19 months of LAT 

514 observing time, and furthermore we find its flux to be consistent with the EGRET measured value, 

515 thus showing no evidence of variability at 7-ray energies over decade timescales as well. The 

516 non- variable 7-ray continuum of the source thus promises the opportunity to observe it over the 

517 extended timescales required for a 5a detection. This is in contrast to the other VHE-detected 

518 quasars, which were detected only during periods where the sources were in a flaring s tate. In this 



519 way 4C +55.17 stands apart from all of the EBL-constraining sources considered in lAbdo et al 



520 (J2010d ). as it holds the greatest potential for providing future constraints. 



521 4. Conclusions 

522 The investigation of the multiwavelength properties of 4C +55.17, including its unusually hard 

523 7-ray spectrum, lack of distinct variability, and CSO-like radio morphology, has highlighted the 

524 exceptional nature of this 7-ray source. For the first time, we have modeled the radio to 7-ray 

525 emission of 4C +55.17 as a young radio source using a dynamic model that is consistent with the 

526 full extent of its observed properties. Furthermore, we find that the prospect of a VHE observation 

527 of 4C +55.17, whose 7-ray spectrum already extends up to the observed energy of 145 GeV, is 

528 within reach of the current generation of Cherenkov telescopes. A detection by such an instrument 

529 would not only add to the present understanding of the source itself, but would also serve to place 



Because the lStecker et all (200a) fast evolution model predicts an increased opacity from the baseline model, our 



conclusions from the baseline test can be applied in both cases. 
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530 a significant upper limit to the level of EBL through combined Fermi and VHE data. Furthermore, 

531 we anticipate that through continued monitoring of 4C +55.17 at high energies with the Fermi 

532 LAT, as well as in the radio through X-rays, the precise classification of 4C +55.17 will become 

533 increasingly more apparent. If, for example, the source continues to remain non-variable in 7-rays 

534 in the years to come, its average fiux versus variability may eventually lie outside the distribution of 

535 Fermi 7-ray emitting blazars altogether, which would make any standard blazar emission scenario 

536 difficult to reconcile. On the other hand, if rapid variability is found in this source, that would 

537 seem to rule out a pure CSO interpretation. Thus we expect that 4C +55.17 will be an important 

538 target for future observations across all wavelengths. 
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561 A. Association of the 145 GeV photon with 4C +55.17 

562 To further investigate the VHE detection of the source, the 145 GeV event was analyzed in 

563 detail using the event displa}0 and found to be a clean 7-ray event, going through more than 

564 half a tracker tower before interacting in the back planes and generating a well-behaved symmetric 

565 shower in the calorimeter. A full Monte Carlo simulation was also run in order to determine the 

566 accuracy of the energy reconstruction. A total of 500, 000 7-rays between the energies 50 and 

567 200 GeV were simulated at an incoming angle 9 and (j) equivalent to that of the measured event. 

568 Data selection cuts were applied on all similar variables, including cuts on the calorimeter raw 

569 energy, best measured energy, reconstructed direction, and event class level. The distribution in 

570 Monte Carlo energy for the remaining events was found to give a ~ Icr error of ±11 GeV. 

571 The probability of the 145 GeV event occurring by random coincidence from background con- 

572 tamination was calculated using the gtsrcprob analysis tool. Probabilities of each event are as- 



573 signe d via standard likelihood analysis to all sources within a provided best-fit model (JMattox et al. 



19961 ). The probability that a photon is produced by a source i is proportional to M\, given by the 



575 formula: 



Mi(e',p',t) 



dedp Si{e,p) R{£',p']e,p,t) , 



(Al) 



SR 



576 where Si{£,p) is the predicted counts density from the source at energy e and position p, and 

577 R{£',p'; e,p, t) is the convolution over the instrument response. In this way, all the surrounding point 

578 sources, the diffuse background, and their corresponding best-fit spectra are taken into account 

579 when assigning probabilities to individual photon events. For the 145 GeV event, the probability of 

580 spurious as sociation with 4 C +5 5.17 was found to be 1.8 x 10~^, agreeing well with an independent 

581 method by lNeronov et al.l (J2010l ). who quote a chance probability by background contamination of 

582 3.1 X 10~^ for the same event. 



B. Calculation of the MAGIC II and VERITAS Differential Flux Sensitivities 



584 Startin g with the integra,l flux sensitivity curves of MAGIC II (JBorla Tridon et al.ll2010l ) and 

585 VERITAS (jPerkins fc Maierl |2009| ) , the differential flux sensitivities can be derived for a given 

586 functional form. In the case of 4C +55.17, we may represent the attenuated VHE spectrum in 

587 general with an exponential cutoff given by the formula: 



dE 



(Bl) 



^http : //glast-ground . slac . Stanford . edu/DataPortalWired/ | 
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588 where Nq, Ec, and F are free parameters of the fitted form of the function. The integral flux above 

589 some minimum energy Eq is thus given by: 

/•oo ^ 

N = No dE E-^ e"^ (B2) 

Jeo 



590 Defining the quantity 



^{E) = / dE E-^ e~^ (B3) 



lEo 
591 the appropriate solution for A^o may be substituted into equation IBll to obtain: 



dN 
'dE 






^{Eo 

-C/O 



(B4) 



592 To construct the differential flux sensitivity curves, we obtained the values F = 2.12 and Ec = 100 

593 GeV by performing a gtlike fit of the > 1.6 GeV data of 4C +55.17 to the exponential cutoff 

594 functional form. For each value N of the integral flux sensitivity, a corresponding differential flux 

595 sensitivity value could thus be obtained via numerical evaluation of equation [ 
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Table 1. Swift/\J\OT Observations of 4C +55.17 



Band 


A [A] 


Fepi [mjy] 


Fcp2 [mJy] 


Fcp'i [mJy] 


V 


5402 


0.331 ±0.061 


0.337 ±0.029 




B 


4329 


0.262 ±0.015 


0.286 ±0.015 




U 


3501 


0.249 ±0.010 


0.251 ±0.011 




UVWl 


2634 


0.175 ±0.007 


0.174 ±0.007 




UVM2 


2231 


0.142 ±0.029 


0.167 ±0.007 




UVW2 


2030 


0.125 ±0.009 


0.127 ±0.005 


0.130 ±0.005 



Note. — The observations were obtained on 2009 March 5 (epl), 
Nov 11 (cp2), and Nov 26 {cp3). 



Table 2. De-absorption of 7-ray flux using different EBL models with fixed t^,^, normalization 



EBL model 



Flux" 



-log(likelihood) 



Finke et al. (2010) 1.83 ± 0.05 2.20 ± 0.06 9.05 ± 0.46 595671.252 

Franceschini et al. (2008) 1.83 ± 0.05 2.21 ± 0.06 9.04 ± 0.46 595671.192 

Gilmore et al. (2009) 1.83 ±0.05 2.21 ±0.06 9.04 ±0.46 595671.133 

Primack et al. (2005) 1.83 ±0.05 2.19 ± 0.06 9.06 ± 0.46 595671.074 

Kneiske (2004) best fit 1.83 ± 0.05 2.18 ± 0.06 9.06 ± 0.46 595671.577 

Kneiske (2004) high UV 1.84 ± 0.05 2.14 ±0.06 9.10 ±0.46 595672.025 

Stecker (2006) baseline 1.85 ± 0.05 2.10 ±0.06 9.14 ±0.46 595678.519 

Stecker (2006) fast evolution 1.85 ± 0.05 2.10 ± 0.07 9.14 ±0.46 595680.170 

Salamon & Stecker (1998) 1.84 ± 0.05 2.15 ± 0.06 9.10 ±0.46 595673.291 



"^Flux above 100 MeV in units of [10^* cm"^ s"^] 
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Fig. 1. — VLBA 5 GHz map (left) featuring the inner parsec-scale radio structure of 4C +55.17, 
reimaged using data from iHelmboldt et al.l (|2007l ). The beam size is 2.0 mas x 1.6 mas (position 
angle = —29.6°), and the contour levels increase by factors of v2 beginning at 1 mJy/beam. The 
resolved morphology has a total angular size of 53 mas (413 pc). The VLA 5 GHz map (right) 
with a 0.4" beam (lo west contour is 2 mJy/ beam increasing by factors of v2) shows the large scale 
radio structure (from iTavecchio et al.l 120071 ). 
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Fig. 2. — Fermi /LAT 19-month 7-ray light curve of 4C +55.17 divided into 7 day bins. All points 
represent > 3a detections and are plotted along with their statistical errors. The dashed horizontal 
line and gray region represent the weighted mean and corresponding error derived from all > 3cj 
detections over the observing period. 
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Fig. 3. — The CSO model of 4C +55.17 versus multiwavelength data, including the new LAT spec- 
trum along with contemporaneou s data with Swift XRT , BAT, and UVOT (black bullets). Archival 
detections (gray) with EGRET (JHartman et al.lll999l ). ROSAT, Chandra, SDSS, 2MASS, 5-year 
integrated WMAP, and historic radio data are also included, as well as archival VLA measurements 
(black triangles) of the inner ~ 400pc radio struct ure (see § 12. 2p . De-absorption of the observed 
Fermi spectral points using the iFinke et al.l (J201Cl ) EBL model was applied in order to properly 
model the intrinsic 7-ray spectrum. Black curves indicate the total non-thermal emission of the 
lobes, with the long-dashed/green representing the contribution from synchrotron self-Compton 
(SSC). Dashed/pink, dash-dot-dotted/gray, and dash-dotted/blue blackbody-type peaks represent 
the dusty torus, starlight, and the UV disk emission components, respectively, along with their 
corresponding inverse-Compton components as required by the model. 
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Fig. 4. — Blazar fit using multi- wavelength data for 4C +55.17. Indicated are the individual con- 
tributions from synchrotron and SSC (long-dashed/green), as well as IC scattering off of the repro- 
cessed UV disk emission from the broad line region (dash-dotted/blue), dusty torus (dashed/pink), 
and host galaxy (dash-dot-dotted/gray); the black curve indicates the total of these components. 
As in Fig.O the dashed/pink, dash-dot-dotted/gray, and dash-dotted/blue blackbody-type peaks 
represent the dusty torus, starlight, and the UV disk emission components, respectively, along with 
their corresponding inverse-Compton components as required by the model. 
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Fig. 5. — The r^^ opacity versus energy for several EBL models ai z = 0.896. The highest- 
energy photon of 145 GeV (rest frame energy = 275 GeV) within the 95% containment radius of 
the 4C +55.17 position is also indicated (vertical dashed line). The horizontal line simply denotes 
r^^ = 1. At the observed energy, attenuation from the EBL is expected even for those models 
which predict low levels of EBL. 



31 



fa 

is 



10 



o 10 



-10 



10 



-11 



10 



-12 



Stacker et al. (2006) fixed t^, normalization 
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Fig. 6. — The de-absorbed spectrum of 4C -1-55.17 (thick da shed hne, gray) along with la error 
bars (thin dashed hnes, red) using the IStecker et al.l (|2006l ) basehne model. Observed spectral 



points without de-absorption, along with the observed spectrum with r^^ normalization left free 
(solid line) and fixed to 1 (dotted line), are plotted for comparison. The de-absorbed spectrum 
shows the non-physical behavior of an unbounded exponential rise up to the observed LAT energy 
of 145 GeV. This trend, which is preferred by 3.9a over a single power law, increases the intrinsic 
spectrum by two orders of magnitude above the observed inverse Compton peak and requires the 
modeling of an additional (and unknown) spectral component (see Figures [3] and H]) . 
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Fig. 7. — The observed LAT spectrum fit to a broken power law with attenuation from 9 different 
EBL models. The spectra are extrapolated beyond the observed energy of 145 GeV and compared 
against the MAGIC II and VERITAS differential flux sensitivity curves for a 50 hour, 5a detection 
of a source characterized by an exponentially decreasing spectrum (see Appendix [B]). For several 
EBL models, the 40 +55.17 spectrum is found to intercept with both the VERITAS and MAGIC II 
sensitivities, making 4C +55.17 a viable candidate for future ground-based VHE observations. 



